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Abstract: We consider SUSY- like missing energy events at hadron colliders and critically 
examine the common assumption that the missing energy is the result of two identical missing 
particles. In order to experimentally test this hypothesis, we generalize the subsystem Mt2 
variable to the case of asymmetric event topologies, where the two SUSY decay chains termi- 
nate in different "children" particles. In this more general approach, the endpoint M T2 ^ max ^ of 
the Mt2 distribution now gives the mass M p (M c (a) , M c (fc) ) of the parent particles as a function 
of two input children masses M c {a) and M c (6) . We propose two methods for an independent 
determination of the individual children masses M c (a) and M c (fe) . First, in the presence of 
upstream transverse momentum Putm the corresponding function M p (Mc a \ Putm) 
is independent of Putm & t precisely the right values of the children masses. Second, the 
previously discussed Mt2 "kink" is now generalized to a "ridge" on the 2-dimensional sur- 
face M p (M c (a) ,M c (b) ). As we show in several examples, quite often there is a special point 
along that ridge which marks the true values of the children masses. Our results allow col- 
lider experiments to probe a multi-component dark matter sector directly and without any 
theoretical prejudice. 
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1. Introduction 

A general expectation in high energy physics today is that physics beyond the standard model 
(BSM) should emerge at the TeV scale in order to stabilize the hierarchy between the Planck 
and electroweak scales. Further evidence in support of this belief is provided by the dark 
matter problem of astro-particle physics, which can be quite naturally solved by postulating 
the existence of a new, weakly-interacting dark matter particle with a mass in the TeV range. 
Such dark matter particles are naturally present in the most popular BSM scenarios such as 
supersymmetry [1], extra dimensions [2-4], little Higgs theory [5,6] etc. They will be produced 
in the upcoming high-energy collisions at the Large Hadron Collider (LHC) at CERN, which 
offers an exciting opportunity to study dark matter in a high-energy lab. Since the dark 
matter particles are weakly interacting, they do not leave any deposits inside the detector 
and can only manifest themselves in the form of missing energy. Recently, there has been a 
lot of theoretical effort directed at testing the dark matter hypothesis at the LHC [7-16] and 
the future International Linear Collider (ILC) [7,9-11,17-20]. Unfortunately, most of these 
studies have been performed in some very model-dependent as well as very complex setup 1 . 
In the literature, a typical collider study of dark matter most often starts with the assumption 
of a specific model with a dark matter candidate (usually supersymmetry with its myriad of 
parameters) and then investigates the model's predictions for the expected rates at the LHC 
in one or several missing energy channels. Rarely, if ever, has the question been posed in 
reverse: what does the observation of a missing energy signal at the LHC tell us about the 
dark matter particle and its properties in a generic and model-independent way [21]. 

1.1 Probing the dark matter sector at colliders 

Naturally, the most pertinent question after the discovery of any BSM missing energy signal 
at the LHC is simply whether the new signal is indeed due to the production of new massive 
particles, or whether it is just an enhancement in the production of SM neutrinos [21]. In 
principle, there are two handles that can be used in addressing this question. In order to 
prove dark matter production, one can measure the mass of the missing particle and show 
that it is different (heavier) from the SM neutrino masses. Alternatively, one can try to 
measure the spin of the missing particle and show that it is different from 1/2 (the spin of 
the neutrino). While there is a large body of recent work on spin measurements in missing 
energy events [22-48], once again very few of those methods are model-independent [45,47]. 
Furthermore, in all considered examples in the literature the spin measurement appears to be 
very difficult. Therefore, in this paper we shall concentrate on the question of measuring the 
mass(es) of the particles responsible for the missing energy. In doing so, we are motivated 
by two reasons. First, previous experience indicates that the mass question will be answered 
long before any spin measurements, and second, many of the spin determination methods 
require prior knowledge of the mass spectrum anyway. 

The difficulty in measuring the mass of the dark matter particle at a hadron collider like 
the Tevatron or the LHC is widely appreciated and has generated a lot of recent activity 

1 Some notable exceptions are the studies in Refs. [7,8]. 
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[49-102]. The main problem can be understood as follows. In a typical BSM dark matter 
scenario, the cosmological longevity of the dark matter particle is ensured by some new 
symmetry 2 under which the SM particles are singlets. At the same time, there are additional 
particles in the spectrum which are charged under the new symmetry. If the lightest one 
among those is electrically and color neutral, it is a potential dark matter candidate, whose 
lifetime is protected by the new symmetry. With any such setup, it is clear that single 
production of dark matter particles at colliders is forbidden by the symmetry Therefore, 
each event has at least two missing particles, whose energies and momenta are unknown. As 
a rule, it is typically impossible to fully reconstruct the kinematics of such events and observe 
the mass of the missing particle directly as an invariant mass peak 3 . Consequently, one has 
to resort to various indirect methods of extracting the mass of the dark matter particle. 

Unfortunately, all existing studies in the literature have explicitly or implicitly made the 
following two assumptions: 

• Single dark matter component. A common assumption throughout the collider phe- 
nomenology literature is that colliders are probing only one dark matter species at a 
time, i.e. that the missing energy signal at colliders is due to the production of one and 
only one type of dark matter particles. Of course, there is no astrophysical evidence 
that the dark matter is made up of a single particle species: it may very well be that 
the dark matter world has a rich structure, just like ours [103]. Consequently, if there 
exist several types of dark matter particles, each contributing some fraction to the total 
relic density, a priori there is no reason why they cannot all be produced in high energy 
collisions. Theoretical models with multiple dark matter candidates have also been 
proposed [104-112]. 

• Identical missing particles in each event. A separate assumption, common to most 
previous studies, is that the two missing particles in each event are identical. This as- 
sumption could in principle be violated as well, even if the single dark matter component 
hypothesis is true. The point is that one of the missing particles in the event may not 
be a dark matter particle, but simply some heavier cousin which decays invisibly. An 
invisibly decaying heavy neutralino (xi ~~ ^ v ^x\ with i > 1) and an invisibly decaying 
sneutrino (y — > vxi) are t wo such examples from supersymmetry. As far as the event 
kinematics is concerned, the mass of the heavier cousin is a relevant parameter and 
approximating it with the mass of the dark matter particle will simply give nonsensical 
results. Another relevant example is provided by models in which the SUSY cascade 
may terminate in any one of several light neutral particles [113]. 

Given our utter ignorance about the structure of the dark matter sector, in this paper we 
set out to develop the necessary formalism for carrying out missing energy studies at hadron 

2 Some popular examples are: imparity in supersymmetry, KK parity in Universal Extra Dimensions, Im- 
parity in Little Higgs models, Z-parity in warped extra dimensions, [/-parity in extended gauge theories, 
etc. 

3 For studies attempting full event reconstruction in long cascade chains, see Refs. [63,75,88,92,101]. 
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colliders in a very general and model-independent way, without relying on any assumptions 
about the nature of the missing particles. In particular, we shall not assume that the two 
missing particles in each event are the same. We shall also allow for the simultaneous pro- 
duction of several dark matter species, or alternatively, for the production of a dark matter 
candidate in association with a heavier, invisibly decaying particle. Under these very general 
circumstances, we shall try to develop a method for measuring the individual masses of all 
relevant particles - the various missing particles which are responsible for the missing energy, 
as well as their parents which were originally produced in the event. 

1.2 Generalizing My 2 to asymmetric event topologies 

In general, by now there is a wide variety of techniques available for mass measurements 
in SUSY- like missing energy events. Such events are characterized by the pair production 
of two new particles, each of which undergoes a sequence of cascade decays ending up in 
a particle which is invisible in the detector. Each technique has its own advantages and 
disadvantages 4 . For our purposes, we chose to revamp the method of the Cambridge Mt2 
variable [50] and adapt it to the more general case of an asymmetric event topology shown 
in Fig. 1. Consider the inclusive production of two identical 5 parents of mass M p as shown 
in Fig. 1. The parent particles may be accompanied by any number of "upstream" objects, 
such as jets from initial state radiation [66,67,79], or visible decay products of even heavier 
(grandparent) particles [85]. The exact origin and nature of the upstream objects will be of 
no particular importance to us, and the only information about them that we shall use will 
be their total transverse momentum Putm- In turn, each parent particle initiates a decay 
chain (shown in red) which produces a certain number of Standard Model (SM) particles 
(shown in gray) and an intermediate "child" particle of mass M C (A) . Throughout this paper 
we shall use the index A to classify various objects as belonging to the upper (A = a) or lower 
(A = b) branch in Fig. 1. The child particle may or may not be a dark matter candidate: 
in general, it may decay further as shown by the dashed lines in Fig. 1. We shall apply the 
"subsystem" Mt2 concept [77, 85] to the subsystem within the blue rectangular frame. The 
SM particles from each branch within the subsystem form a composite particle of known 6 
transverse momentum and invariant mass m<\\. Since the children masses and 
are a priori unknown, the subsystem Mt2 will be defined in terms of two "test" masses 
Mc°^ and In Fig. 1, q^ are the trial transverse momenta of the two children. The 

individual momenta q^ are also a priori unknown, but they are constrained by transverse 
momentum conservation: 

q^ + = Q tot = -(p^ + + p UTM ). (1.1) 

Given this very general setup, in Section 3 we shall consider a generalization 7 of the usual 
4 For a comparative review of the three main techniques, see [85]. 

5 In principle, the assumption of identical parents can also be relaxed, by a suitable generalization of the 
Mt2 variable, in which the mass ratio of the two parents is treated as an additional input parameter [96]. 

6 We assume that there are no neutrinos among the SM decay products in each branch. 

7 The possibility of applying the Mt2 variable to an event topology with different children was previously 
mentioned in Refs. [95,96]. 
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Figure 1: The generic event topology under consideration in this paper. We consider the inclusive 
pair-production of two "parent" particles with identical masses M p . The parents may be accompanied 
by "upstream" objects, e.g. jets from initial state radiation, visible decay products of even heavier 
particles, etc. The transverse momentum of all upstream objects is measured and denoted by Pjjtm- 
In turn, each parent particle initiates a decay chain (shown in red) which produces a certain number 



of SM particles (shown in gray) and an intermediate "child" particle of mass Mj , where X 



(X = b) for the branch above (below). In general, the child particle does not have to be the dark matter 
candidate, and may decay further as shown by the dashed lines. The Mt2 variable is defined for the 
subsystem inside the blue box and is defined in terms of two arbitrary children "test" masses M c (a) ana 



M, 



(6) 



The SM particles from each branch form a composite particle of transverse momentum pj, 



(A) 



and invariant mass m(X) > correspondingly. The trial transverse momenta q>p of the children obey 
the transverse momentum conservation relation shown inside the green box. In general, the number 
n^ x \ as well as the type of SM decay products in each branch do not have to be the same. 



Mt2 variable which can apply to the asymmetric event topology of Fig. 1. There will be two 
different aspects of the asymmetry: 



• First and foremost, we shall avoid the common assumption that the two children have 
the same mass. This will be important for two reasons. On the one hand, it will allow 
us to study events in which there are indeed two different types of missing particles. 
We shall give several such examples in the subsequent sections. More importantly, the 
endpoint of the asymmetric Mt2 variable will allow us to measure the two children 
masses separately. Therefore, even when the events contain identical missing particles, 
as is usually assumed throughout the literature, one would be able to establish this fact 
experimentally from the data, instead of relying on an ad hoc theoretical assumption. 
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• As can be seen from Fig. 1, in general, the number as well as the types of SM decay 
products in each branch may be different as well. Once we allow for the children to be 
different, and given the fact that we start from identical parents, the two branches of 
the subsystem will naturally involve different sets of SM particles. 

In what follows, when referring to the more general Mt2 variable defined in Section 3, we 
shall interchangeably use the terms "asymmetric" or "generalized" Mt2- In contrast, we 
shall use the term "symmetric" when referring to the more conventional Mt2 definition with 
identical children. 

The traditional Mt2 approach assumes that the children have a common test mass M c = 
= and then proceeds to find one functional relation between the true child mass 
M c and the true parent mass M p as follows [50]. Construct several Mt2 distributions for 
different input values of the test children mass M c and then read off their upper kinematic 
endpoints M T2 t max \(M c ). These endpoint measurements are then interpreted as an output 
parent mass M p , which is a function of the input test mass M c : 

M p (M c ) = M T2(max) (M c ) . (1.2) 

The importance of this functional relation is that it is automatically satisfied for the true 
values M p and M c of the parent and child masses: 

M p = M T2(max) (M c ). (1.3) 

In other words, if we could somehow guess the correct value M c of the child mass, the function 
(1.2) will provide the correct value M p of the parent mass. However, since the true child mass 
M c is a priori unknown, the individual masses M p and M c still remain undetermined and 
must be extracted by some other means. 

At this point, it may seem that by considering the asymmetric Mt2 variable with non- 
identical children particles, we have regressed to some extent. Indeed, we are introducing an 
additional degree of freedom in eq. (1.2), which now reads 

M p (MW,MW) = M T2{max) (M^,M^) . (1.4) 

The standard M?2 endpoint method will still allow us to find the parent mass M p , but now 
it is a function of two input parameters and which are completely unknown. Of 

course, if one knew the correct values of the two children masses M c (a) and M c (b) entering 
eq. (1.4), the true parent mass M p will be given in a manner analogous to eq. (1.3): 

M p = M T2{max) {M^,M^). (1.5) 

Our main result in this paper is that in spite of the apparent remaining arbitrariness 
in eq. (1.4), one can nevertheless uniquely determine all three masses M p , and 
just by studying the behavior of the measured function M p {Mc a \M^). More importantly, 
this determination can actually be done in two different ways! Our first method is simply a 
generalization of the observation made in Refs. [65-68, 85] that under certain circumstances 
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(varying fntX) or nonvanishing upstream momentum Pjjtm), the function (1.2) develops a 
"kink" precisely at the correct value M c of the child mass: 

' dM p {M c ) \ ( dM p (M c ) \ J = 0, if M c + M c , 



\ dM c I - V dM c ~ / 0, if M c = M c . 

In other words, the function (1.2) is continuous, but not differentiable at the point M c = M c . 
In the asymmetric Mt2 case, we find that the function (1.4) is similarly non-differentiable 
at a set of points {(Mc a \ M^)}, so that the kink of eq. (1.6) is generalized to a "ridge" 
on the 2-dimensional hypersurface defined by (1.4) in the three-dimensional parameter space 
of {M c (a) , M c (6) , M p } . 8 Interestingly enough, the ridge often (albeit not always) exhibits a 
special point which marks the exact location of the true values (M^ , ) . 

Our second method for determining the two children masses M c (a) and M c (fe) is even more 
general and is applicable under any circumstances. The main starting point is that just like 
the endpoint of the symmetric Mt2, the endpoint of the asymmetric Mt2 also depends on 
the value of the upstream transverse momentum Putm, so that eq. (1.4) is more properly 
written as 

M P (M^ , , P UTM ) = M T2(max) (M^ , Mf) , P UTM ) . (1.7) 

Now we can explore the Pjjtm dependence in (1.7) and note that it is absent for precisely 
the right values of and : 



dM T2(max) {Mi a \M^ b \p UTM ) 



dPi 



UTM 



While this property has been known, it was rarely used in the case of the symmetric Mt2 , since 
it offers redundant information: once the correct child mass M c is found through the Mt2 
kink (1.6), the parent mass M p is given by (1.2) and there are no remaining unknowns, thus 
there is no need to further investigate the Pjjtm dependence. In the case of the asymmetric 
Mt2, however, we start with one additional unknown parameter, which cannot always be 
determined from the "ridge" information alone. Therefore, in order to pin down the complete 
spectrum, we are forced to make use of (1.8). The nice feature of the Putm method is that 
it always allows us to determine both children masses and , without relying on the 
"ridge" information at all. In this sense, our two methods are complementary and each can 
be used to cross-check the results obtained by the other. 

The paper is organized as follows 9 . In Sec. 2 we begin with a review of the conventional 
symmetric Mt2 variable and its properties. Then in Sec. 3 we introduce the asymmetric Mt2 
variable and highlight its properties which are relevant for our mass measurements. We also 
discuss some experimental subtleties in the construction of the asymmetric Mt2 distribution, 
which are not present in the case of the symmetric Mt2- Sections 4, 5.1 and 5.2 present some 



8 Ref. [96] studied the orthogonal scenario of different parents (Mp"' 7^ M^) and identical children (jVfj a ' = 
Mc b ^) and found a similar non-differentiable feature, called a "crease", on the corresponding two-dimensional 
hypersurface within the three-dimensional parameter space 

9 Readers who are unfamiliar with the Mt2 concept may benefit from consulting Refs. [54,68,85,96] first. 
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simple examples of asymmetric event topologies. Finally, Sec. 6 summarizes our main results 
and outlines some possible directions for future work. Appendix A revisits the examples of 
Section 4 in the case of Putm — > co, which can be handled by purely analytical means [96]. 



2. The conventional symmetric M T2 
2.1 Definition 

We begin our discussion by revisiting the conventional definition of the symmetric Mt2 
variable with identical daughters, following the general notation introduced in Fig. 1. Let us 
consider the inclusive production of two parent particles with common mass M p . Each parent 
initiates a decay chain producing a certain number of SM particles. In this section we 
assume that the two chains terminate in children particles of the same mass: M c (a) = M c (fe) = 
M c . (From Section 3 on we shall remove this assumption.) In most applications of Mt2 in 
the literature, the children particles are identified with the very last particles in the decay 
chains, i.e. the dark matter candidates. However, the symmetric M?2 can also be usefully 
applied to a subsystem of the original event topology, where the children are some other pair 
of (identical) particles appearing further up the decay chain [77,85]. The Mt2 variable is 
defined in terms of the measured invariant mass m^) an d transverse momentum of the 
visible particles on each side (see Fig. 1). With the assumption of identical children, the 
transverse mass of each parent is 



(p T (A) ; <f T (A) ; m (A) ; M c ) = Jmf x) + Mf + 2 ( e WgW -p T W • <f T (A) ) , (2.1) 



where M c is the common test mass for the children, which is an input to the Mt2 calculation, 
while is the unknown transverse momentum of the child particle in the A-th chain. 

In eq. (2.1) we have also introduced shorthand notation for the transverse energy of the 
composite particle made from the visible SM particles in the A-th chain 



e^^+^'rf 1 (2.2) 
and for the transverse energy of the corresponding child particle in the A-th chain 

e™=y/M*+q T W -q T W . (2.3) 

Then the event-by-event symmetric M?2 variable is defined through a minimization pro- 
cedure over all possible partitions of the two children momenta q^ [50] 



Mt2 (pt a) <Pt 6) ' m (a) M c , P UT m) 



mm 

_ (a) . ^ (6) p, 

<? T V ' +q T -Qtot 



max 



[M^(p^; q T ^; m (a) ; M c ) , (p^ ; m (b) ; M c ) } 



, (2-4) 



consistent with the momentum conservation constraint (1.1) in the transverse plane. 
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2.2 Computation 

The standard definition (2.4) of the Mt2 variable is sufficient to compute the value of Mt2 
numerically, given a set of input values for its arguments. The right-hand side of eq. (2.4) 
represents a simple minimization problem in two variables, which can be easily handled by a 
computer. In fact, there are publicly available computer codes for computing Mt2 [114,115]. 
The public codes have even been optimized for speed [84] and give results consistent with each 
other (as well as with our own code) 10 . Nevertheless, it is useful to have an analytical formula 
for calculating the event-by-event Mt2 for several reasons. First, an analytical formula is 
extremely valuable when it comes to understanding the properties and behavior of complex 
mathematical functions like (2.4). Second, computing Mt2 from a formula will be faster than 
any numerical scanning algorithm. The computing speed becomes an issue especially when 
one considers variations of Mt2 like MT2gen, where in addition one needs to scan over all 
possible partitions of the visible objects into two decay chains [64]. Therefore, in this paper 
we shall pay special attention to the availability of analytical formulas and we shall quote 
such formulas whenever they are available. 

In the symmetric case with identical children, an analytical formula for the event-by- 
event Mt2 exists only in the special case Pjjtm = 0. It was derived in [64] and we provide 
it here for completeness. (In the next section we shall present its generalization for the 
asymmetric case of different children.) The symmetric Mt2 is known to have two types of 
solutions: "balanced" and "unbalanced" [54,64]. The balanced solution is achieved when the 
minimization procedure in eq. (2.4) selects a momentum configuration for in which the 
transverse masses of the two parents are the same: = In that case, typically 

neither nor M^p is at its global (unconstrained) minimum. In what follows, we shall 

use a superscript B to refer to such balanced-type solutions. The formula for the balanced 
solution My 2 of the symmetric Mt2 variable is given by [64, 68] 



^T2(PT a) 'PT b) ' m (a).™-(6); Mc) 

(2.5) 

where At is a convenient shorthand notation introduced in [68] 

A T = e^e^ +p^ a) -p^ b) (2.6) 

and was already defined in eq. (2.2). 

On the other hand, unbalanced solutions arise when one of the two parent transverse 
masses (M^ a) orM?, as the case may be) is at its global (unconstrained) minimum. Denoting 
the two unbalanced solutions with a superscript UX, we have [54] 

Mt2 (p T (a) ,P T (b) ; "i(a) ,m ( 6) ; M c ) = m {a) + M C , (2.7) 

Mt2 {Pr a) >PT b) ; m (a) ,m (6) ; M c ) = m (6) + M c . (2.8) 

10 Unfortunately, the assumption of identical children is hardwired in the public codes and they cannot be 
used to calculate the asymmetric Mt2 introduced below in Section 3 without additional hacking. We shall 
return to this point in Section 3. 



= M^+A T + 



\ 



AM* 



2 At 



m 



(a) 



m 



(b), 



-Arjn 777/ / \ 7Tl/i\ 

T (a) (b) 
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Given the three possible options for Mt2, eqs. (2.5), (2.7) and (2.8), it remains to specify 
which one actually takes place for a given set of values for p^ a \ p^ b \ fn< a ), m (b), M c and 
Putm = in the event 11 . The balanced solution (2.5) applies when the following two 
conditions are simultaneously satisfied: 



(p T W; ? T W =-* T $+ Qtof, m (b y, Mo) > M^ a) (p T M; g T (a) =g T %\; m (o) , M c ) = m (a) + M c , (2.9) 

M^ a) (p T W; ? T W =-?^+ Qtof, m (a) ; M c ) > M^(p T (6) ; q T W =q T %; "W = m (*>) + M c , (2-10) 

where 

(A = a,6), (2.11) 



(a) _ M, 

?T(0) " m 



c -> (A) 
(A) 



(A) 



gives the global (unconstrained) minimum of the corresponding parent transverse mass M 
The unbalanced solution M^g applies when the condition (2.9) is false and condition (2.10) is 
true, while the unbalanced solution applies when the condition (2.9) is true and condition 
(2.10) is false. It is easy to see that conditions (2.9) and (2.10) cannot be simultaneously 
violated, so these three cases exhaust all possibilities. 

2.3 Properties 

Given its definition (2.4), one can readily form and study the differential Mt2 distribution. 
Although its shape in general does carry some information about the underlying process, 
it has become customary to focus on the upper endpoint M T2 ( ma x), which is simply the 
maximum value of Mt2 found in the event sample: 



M 



T2(max) (M c , Putm) = max 

all events 



MT2 (p T (<2) ,PT b) ; "l( o) ,m( 6 ); A*c) 



(2.12) 



Notice that in the process of maximizing over all events, the dependence on , p^ , 
and mn,) disappears, and Mx2{max) depends only on two input parameters: M c and Putm, the 
latter entering through Qtot in the momentum conservation constraint (1.1). The measured 
function (2.12) is the starting point of any My2-based mass determination analysis. We shall 
now review its three basic properties which make it suitable for such studies [100]. 

2.3.1 Property I: Knowledge of M p as a function of M c 

This property was already identified in the original papers and served as the main motivation 
for introducing the Mt2 variable in the first place [50,54]. Mathematically it can be expressed 
as 



M p (M c ,P UTM ) = M T2(max) (M c ,P UTM ). 



(2.13) 



This is the same as eq. (1.2), but now we have been careful to include the explicit dependence 
on Putm, which will be important in our subsequent discussion. As indicated in eq. (2.13), 
the function M p (M c , Putm) can be experimentally measured from the Mt2 endpoint (2.12). 



L Recall that (2.5) only applies for Putm = 0. 
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The crucial point now is that the relation (2.13) is satisfied by the true values M p and M c of 
the parent and child mass, correspondingly: 

M P = M T2{max) (M c , P UTM ) - (2.14) 

Notice that eq. (2.14) holds for any value of Putm, so in practical applications of this method 
one could choose the most populated Pjjtm bin to reduce the statistical error. On the other 
hand, since a priori we do not know the true mass M c of the missing particle, eq. (2.14) 
gives only one relation between the masses of the mother and the child. This is illustrated in 
Fig. 2(a), where we consider the simple example of direct slepton pair production 12 , where 
each slepton (£) decays to the lightest neutralino (x?) by emitting a single lepton I: I — > £+Xi- 
Here the slepton is the parent and the neutralino is the child. Their masses were chosen to be 
M p = 300 GeV and M c = 100 GeV, correspondingly, as indicated with the black dotted lines 
in Fig. 2(a). In this example, the upstream transverse momentum Pjjtm is provided by jets 
from initial state radiation. In Fig. 2(a) we plot the function (2.13) versus M c , for several fixed 
values of Putm- The green solid line represents the case of no upstream momentum Pjjtm = 
0. In agreement with eq. (2.14), this line passes through the point (M c , M p ) corresponding to 
the true values of the mass parameters. Notice that the property (2.14) continues to hold for 
other values of Pjjtm- Fig. 2(a) shows three more cases: Pjjtm = 500 GeV (dotdashed black 
line), Pjjtm = 1 TeV (dashed red line) and Pjjtm = 2 TeV (dotted blue line). All those 
curves still pass through the point (M c , M p ) with the correct values of the masses, illustrating 
the robustness of the property (2.14) with respect to variations in Pjjtm- 

2.3.2 Property II: Kink in M T2 (max) a * the true M c 

The second important property of the Mt2 variable was identified rather recently [65-68,85]. 
Interestingly, the Mt2 endpoint Mf2( mm ), when considered as a function of the unknown 
input test mass M c , often develops a kink (1.6) at precisely the correct value M c = M c of the 
child mass. The appearance of the kink is a rather general phenomenon and occurs under 
various circumstances. It was originally noticed in event topologies with composite visible 
particles, whose invariant mass rnr\\ is a variable parameter [65,68]. Later it was realised 
that a kink also occurs in the presence of non-zero upstream momentum Putm [66,67,85], as 
in the example of Fig. 2(a), where Pjjtm arises due to initial state radiation. As can be seen 
in the figure, the kink is absent for Pjjtm = 0, but as soon as there is some non-vanishing 
Pjjtm , the kink becomes readily apparent. As expected, the kink location (marked by the 
vertical dotted line) is at the true child mass (M c = 100 GeV), where the corresponding value 
of M T2 (max) (marked by the horizontal dotted line) is at the true parent mass (M p = 300 
GeV). Fig. 2(a) also demonstrates that with the increase in Putm, the kink becomes more 
pronounced, thus the most favorable situations for the observation of the kink are cases with 
large Putm, e.g. when the upstream momentum is due to the decays of heavier (grandparent) 
particles [85]. 

12 The corresponding event topology is shown in Fig. 3(a) below with M^"' = = M c . 
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Figure 2: Plots of (a) the Mt2 endpoint M T2 (max){M c , Putm) defined in eq. (2.12), and (b) the 
function AMT2(max)(M c , Putm) defined in (2.16) as a function of the test child mass M c , for several 
fixed values of Putm-' Putm = GeV (solid, green), Putm = 500 GeV (dot-dashed, black), Putm 
1 TeV (dashed, red), and Putm = 2 TeV (dotted, blue). The process under consideration is pair 
production of sleptons of mass M p = 300 Ge V, which decay directly to the lightest neutralino of 
mass M c = 100 GeV. 



In Sec. 3.3 we shall see how the kink feature (1.6) of the symmetric Mt2 endpoint 
M p (M c ) defined by eq. (1.2) is generalized to a "ridge" feature on the asymmetric M?2 
endpoint M p (M c (a) , M c (b) ) defined in (1.4). 

2.3.3 Property III: Pjjtm invariance of M T2 ( max ) at the true M c 

This property is the one which has been least emphasized in the literature. Notice that the 
Mt2 endpoint function (2.13) in general depends on the value of Pjjtm- However, the first 
property (2.14) implies that the Pjjtm dependence disappears at the correct value M c of the 
child mass: 

dM T2{max) (M c ,P UTM ) 



= 0. (2.15) 

M C =M C 



dPjJTM 

In order to quantify this feature, let us define the function 

AM T2(max) (M c , P UTM ) = M T2{max) (M c , P UTM ) - M T2(max) (M c , 0), (2.16) 

which measures the shift of the Mt2 endpoint due to variations in Pjjtm- The function 
/S.M T2 (max) (M c , Pjjtm) can be measured experimentally: the first term on the right-hand 
side of (2.16) is simply the Mt2 endpoint observed in a subsample of events with a given 
(preferably the most common) value of Putm, while the second term on the right-hand side of 
(2.16) contains the endpoint Mj^^ of the 1-dimensional Mt2 ± variable introduced in [100]: 

M T2(max) (M c ,0) = M^f(M c ). (2.17) 



- 12 - 



Given the definition (2.16), the third property (2.15) can be rewritten as 

AM T2(max) (M c ,P UTM )>0, (2.18) 

where the equality holds only for M c = M c : 

AM T2(max) (M c , P UTM ) = 0, VPutm. (2.19) 

Eqs. (2.18) and (2.19) provide an alternative way to determine the true child mass M c : simply 
find the value of M c which minimizes the function AMr2( max )(M C) Putm)- This procedure 
is illustrated in Fig. 2(b), where we revisit the slepton pair production example of Fig. 2(a) 
and plot the function AM T2 / max \(M c , Putm) defined in (2.16) versus the test mass M c , for 
the same set of (fixed) values of Putm- Clearly, the zero of the function (2.16) occurs at 
the true child mass M c = M c = 100 GeV, in agreement with eq. (2.19). In our studies of 
the asymmetric Mt2 case in the next sections, we shall find that the third property (2.19) 
is extremely important, since it will always allow us the complete determination of the mass 
spectrum, including both children masses M c (a) andM c (6) . 



3. The generalized asymmetric M 



T2 



After this short review of the basic properties of the conventional symmetric Mt2 variable 
(2.4), we now turn our attention to the less trivial case of mJ°^ / Following the logic 

of Sec. 2, in Sec. 3.1 we first introduce the asymmetric Mt2 variable and then in Sees. 3.2 
and 3.3 we discuss its computation and mathematical properties, correspondingly. 

3.1 Definition 

The generalization of the usual definition (2.4) to the asymmetric case of M c (a) / M c (6) is 
straightforward [96]. We continue to follow the conventions and notation of Fig. 1, but now 
we simply avoid the assumption that the children masses are equal, and we let each one be 
an independent input parameter ■ Without loss of generality, in what follows we assume 
M { c ] > M c (a) . The transverse mass of each parent (2.1) is now a function of the corresponding 
child mass Mj^: 

M T X) (p^ X) ; q^; m (A) ; AfW) = ^ + (m c (A) ) 2 + 2 (eWgW - • q^) , (3.1) 

where the transverse energy of the composite SM particle on the A-th side of the event 
was already defined in (2.2), while the transverse energy e^ A ^ of the child is now generalized 
from (2.3) to 



Mr) +qr-QT- (3-2) 
The event-by-event asymmetric Mt2 variable is defined in analogy to (2.4) and is given by [96] 

M T2 (p T (a \p T (b) ;m (a) ,m ib yM a) M b \p UTM ) = 



mm 

_(o),_(6) n 

1 T y ' -Qtot 



max 



{M^M; m (a) ; Af<«>) , M ? (f™ ; ; m (b y, M™)} 



, (3.3) 
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which is now a function of two input test children masses M c (a) and M c (fe) . In the special 
case of mJ°^ = M^ = M c , the asymmetric Mt2 variable defined in (3.3) reduces to the 
conventional symmetric Mf2 variable (2.4). 

3.2 Computation 

In this subsection we generalize the discussion in Section 2.2 and present an analytical formula 
for computing the event-by-event asymmetric Mt2 variable (3.3). Just like the formula (2.5) 
for the symmetric case, our formula will hold only in the special case of Pjjtm = 0. As 
before, the asymmetric Mt2 variable has two types of solutions - balanced and unbalanced. 
The balanced solution occurs when the following two conditions are simultaneously satisfied 
(compare to the analogous conditions (2.9) and (2.10) for the symmetric case) 

M^ ) (? T W;f T W=-f T W+ Qtof,m (b yM b) ) > M^\p^^ a) =q T %,m {a) ,M^) = ro (o) + M c (a) , (3.4) 

M T a) (p T ^;g T ^=-q T %+ Q toi ;m (a) ;M C W ) > (p T (b) ;<? T (b) =<f T %;m (b) M b) ) = m {b) + , (3.5) 



where, in analogy to (2.11) 



-> (A) M c ^( A ) . . . . 

1t)o) = Pt > (A = 0,6), 3.6 



is the test child momentum at the global unconstrained minimum of M^\ The balanced 
solution for Mt2 is now given by 



M&faMtfV; m (o) ,m (6) ; Af<"> M b) ) 



= Ml + A T+ I ® V )mI 
yiA T -m (a) -m (b) J 



± 



4M? { 2M 2 



2 



1 + — ...2 + ...2 +1— —2 _ m 2 1 X \/4- m L™L. ( 3 - 7 ) 

(a) m (b) , 



\ 2A T -m; a) -m; 6) T ^ T -mL-m 2 . / ' V ">>"\M- 



where At was defined in (2.6). For convenience, in (3.7) we have introduced two alternative 
mass parameters 



1 f / ,~r(h)\ 2 



Ml = l\(M^Y + (Mto) 2 }, (3.8) 



^ = l[(M^) 2 - (Mto) 2 } , (3.9) 

in place of the original trial masses mJ°^ and M^\ The new parameters M+ and M_ 
are simply a different parametrization of the two degrees of freedom corresponding to the 
unknown child masses M^ and M^ entering the definition of the asymmetric Mt2- The 
parameters M + and M_ allow us to write formula (3.7) in a more compact form. More 
importantly, they also allow to make easy contact with the known results from Section 2 by 
taking the symmetric limit M^ = M^ = M c as 

M + -»• M c , M_ -»• 0. (3.10) 
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It is easy to see that in the symmetric limit (3.10) our balanced solution (3.7) for the asym- 
metric Mt2 reduces to the known result (2.5) for the symmetric Mt2- 

An interesting feature of the asymmetric balanced solution is the appearance of a ± sign 
on the second line of (3.7). In principle, this sign ambiguity is present in the symmetric case 
as well, but there the minus sign always turns out to be unphysical and the sign issue does 
not arise [64]. However, in the asymmetric case, both signs can be physical sometimes and 
one must make the proper sign choice in eq. (3.7) as follows. For the given set of test masses 
(Mj a ^ , ) , calculate the transverse center-of-mass energy 

± _ 2 (eW - e('>)Ml , (e W + e^)A T - (e^m^ + e^m\ h) ) 

2A T - m \ a) - m\ b) I j% _ 



AMI ( 2M 2 \ 

1 + ~A ^2 + -2 +1-A „2~ m 2 > (3-H) 

(a) m (b) ) 



corresponding to each sign choice in eq. (3.7), and compare the result to the minimum allowed 
value of V§ T 



^ T(mm) = e W + e <»> + J C? M + (m™ + M®)* • (3.12) 

The minus sign in eq. (3.7) takes precedence and applies whenever it is physical, i.e. whenever 
s/i T > V$ T (min)- I n the remaining cases when \f§ T < V§ T (min) an d the minus sign is 
unphysical, the plus sign in eq. (3.7) applies. 

If one of the conditions (3.4), (3.5) is not satisfied, the asymmetric Mt2 is given by an 
unbalanced solution, in analogy to (2.7) and (2.8): 

M^(pT (a W 6) ; m (a) ,m (b y, M») = m (a) + , (3.13) 

M$(pM#W; m (o)l m (6) ; M^,M^) = m (b) + . (3.14) 

The unbalanced solution Mj£ of eq. (3.13) applies when the condition (3.4) is false and 
condition (3.5) is true, while the unbalanced solution of eq. (3.14) applies when the 

condition (3.4) is true and condition (3.5) is false. 

Eqs. (3.7), (3.13) and (3.14) represent one of our main results. They generalize the 
analytical results of Refs. [64, 68] and allow the direct computation of the asymmetric Mt2 
variable without the need for scanning and numerical minimizations. This is an important 
benefit, since the existing public codes for Mt2 [114,115] only apply in the symmetric case 
M c (a) = M c (b) . 

3.3 Properties 

All three properties of the symmetric Mt2 discussed in Section 2.3 readily generalize to the 
asymmetric case. 
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3.3.1 Property I: Knowledge of M p as a function of and 

In the asymmetric case, the endpoint M^max) of the Mt2 distribution still gives the mass 
of the parent, only this time it is a function of two input test masses for the children: 

M p (M^ , , P UTM ) = M T2(max) (M C W , , P UTM ) • (3.15) 

The important property is that this relation is satisfied by the true values of the children and 
parent masses: 

M p = M T2{max) (M^,M^,P UTM ) . (3.16) 

Thus the true parent mass M p will be known once we determine the two children masses 
M c (a) and M c (b) . 

3.3.2 Property II: Ridge in M T2(jriax) through the true M c (a) and M c (b) 

In the symmetric M^ 2 case, the endpoint function (2.13) is not continuously differentiable 
and has a "kink" at the true child mass M c = M c . In the asymmetric Mt2 case, the endpoint 
function (3.15) is similarly non-differentiable at a set of points 

{{m^(9\M^{9))} (3.17) 

parametrized by a single continuous parameter 9. The gradient of the endpoint function 
(3.15) suffers a discontinuity as we cross the curve defined by (3.17). Since (3.15) represents 
a hypersurface in the three-dimensional parameter space of {M c (a) ,M c (b) ,Mp}, the gradient 
discontinuity will appear as a "ridge" (sometimes also referred to as a "crease" [96]) on our 
three-dimensional plots below. The important property of the ridge is that it passes through 
the correct values for the children masses, even when they are different: 

AfW = M c (a) (# ), (3.18) 
= M c (fe) (#o), (3.19) 

for some 9q. Thus the ridge information provides a relation among the two children masses 
and leaves us with just a single unknown degree of freedom — the parameter 9 in eq. (3.17). 

Interestingly, the shape of the ridge provides a quick test whether the two missing par- 
ticles are identical or not 13 . If the shape of the ridge in the (mJ°\ Mc^) plane is symmetric 
with respect to the interchange mJ°^ <H> i.e. under a mirror reflection with respect to 

the 45° line M c (a) = M c (b) , then the two missing particles are the same. Conversely, when the 
shape of the ridge is not symmetric under ^ the missing particles are in general 

expected to have different masses. 



3 To be more precise, the ridge shape tests whether the two missing particles have the same mass or not. 
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3.3.3 Property III: Pjjtm invariance of M T2 ( max ) at the true and 

The third Mt2 property, which was discussed in Section 2.3.3, is readily generalized to the 
asymmetric case as well. Note that eq. (3.16) implies that the Pjjtm dependence of the 
asymmetric Mt2 endpoint (3.15) disappears at the true values of the children masses: 



dM T2{max) {M^ a \M^\p UTM ) 



dPi 



UTM 



M^=Mi a \M^=M^ ' (3 ' 20) 



This equation is the asymmetric analogue of eq. (2.15). Proceeding as in Sec. 2.3.3, let us 
define the function 

AM T2{max) (M^,M( b \P UTM ) = M T2{max) (M^,M^,P UTM ) - M T2{max) (M^,M( b \0), 

(3.21) 

which quantifies the shift of the asymmetric Mt2 endpoint (3.15) in the presence of Putm- 
By definition, 

AM T2(m(tt) (M C W , , Putm) > 0, (3.22) 
with equality being achieved only for the correct values of the children masses: 

AM T2(mflx) (M« , M C W , P UTM ) = 0, V P UTM . (3.23) 

The last equation reveals the power of the Putm invariance method. Unlike the kink method 
discussed in Sec. 3.3.2, which was only able to find a relation between the two children 
masses M c (a) and M c (b) , the Putm invariance implied by eq. (3.23) allows us to determine 
each individual children mass, without any theoretical assumptions, and even in the case 
when the two children masses happen to be different (M c (a) / M c (6) ^ 



3.4 Examples 

In the next two sections we shall illustrate the three properties discussed so far in Section 3.3 
with some concrete examples. Instead of the most general event topology depicted Fig. 1, 
here we limit ourselves to the three simple examples shown in Fig. 3. 

The simplest possible case is when = 1, i.e. when each cascade decay contains a single 
SM particle, as in Fig. 3(a). In this example, rri(x) is constant. For simplicity, we shall take 
m (\) ~ 0) which is the case for a lepton or a light flavor jet. If the SM particle is a Z-boson 
or a top quark, its mass cannot be neglected, and one must keep the proper value of m^y 
This, however, is only a technical detail, which does not affect our main conclusions below. 
In spite of its simplicity, the topology of Fig. 3(a) is actually the most challenging case, due 
to the limited number of available measurements [85]. In order to be able to determine all 
individual masses in that case, one must consider events with upstream momentum Putm, 
as illustrated in Fig. 3(a). This is not a particularly restrictive assumption, since there is 
always a certain amount of Putm in the event (at the very least, from initial state radiation). 
In Section 4 the topology of Fig. 3(a) will be extensively studied - first for the asymmetric 
case of / in Sec. 4.1, and then for the symmetric case of = in Sec. 4.2. 
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(a) m (6) =0 (b) m (6) (c) m (6) 

Figure 3: TTie i/iree different event-topologies under consideration in this paper. In each case, two 
parents with mass M p are produced onshell and decay into two daughters of (generally different) masses 
Mc and Mc ■ Case (a), which is the subject of Section 4, has a single massless visible SM particle 
in each leg and some arbitrary upstream transverse momentum Pjjtm- In the remaining two cases (b) 
and (c), which are discussed in Section 5, there are two massless visible particles in each leg, which 
form a composite visible particle with varying invariant mass m(\) . The intermediate particle of mass 
M^ X) is (b) heavy and off-shell (m[ X) > M p ), or (c) on-shell (M p > M^ X] > m[ X) ). For simplicity, 
we do not consider any upstream momentum in cases (b) and (c). 



Another simple situation arises when there are two massless visible SM particles in each 
leg, as illustrated in Figs. 3(b) and 3(c). In either case, the invariant mass ?Ti(A) is n °t constant 

q^ x , where m™" 



any more, but varies within a certain range m7£™ < < m™ x , where m m ' 



0, while 

the value of mm 1 depends on the mass of the corresponding intermediate particle. In 

Fig. 3(b) we assume > M p , so that the intermediate particle is off-shell and 



(A) 



(3.24) 



The "off-shell" case of Fig. 3(b) will be discussed in Sec. 5.1. 

In contrast, in Fig. 3(c) we take M p > m\ X) > M C (A) , in which case the intermediate 
particle is on-shell and the range for m^x) is now limited from above by 



m 



(A) 



M n 



M (A) V 



M, 



(A)' 



(A) 



(3.25) 



We shall discuss the "on-shell" case of Fig. 3(c) in Sec. 5.2. 

In the event topologies of Figs. 3(b) and 3(c), the mass is varying and the ridge 
of eq. (3.17) will appear even if there were no upstream transverse momentum in the event. 
Therefore, in our discussion of Figs. 3(b) and 3(c) in Sec. 5 below we shall assume Pjjtm — 
for simplicity. The presence of non-zero Pjjtm will only additionally enhance the ridge 
feature. 

3.5 Combinatorial issues 

Before going on to the actual examples in the next two sections, we need to discuss one minor 
complication, which is unique to the asymmetric M?2 variable and was not present in the 
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case of the symmetric Mt2 variable. The question is, how does one associate the visible decay 
products observed in the detector with a particular decay chain A = a or A = b. This is the 
usual combinatorics problem, which now has two different aspects: 

• The first issue is also present in the symmetric case, where one has to decide how to 
partition the SM particles observed in the detector into two disjoint sets, one for each 
cascade. In the traditional approach, where the children particles are assumed to be 
identical, the two sets are indistinguishable and it does not matter which one is first 
and which one is second. This particular aspect of the combinatorial problem will also 
be present in the asymmetric case. 

• In the asymmetric case, however, there is an additional aspect to the combinatorial 
problem: now the two cascades are distinguishable (by the masses of the child particles), 
so even if we correctly divide the visible objects into the proper subsets, we still do not 
know which subset goes together with and thus gets a label A = a, and which goes 
together with and gets labelled by A = b. This leads to an additional combinatorial 
factor of 2 which is absent in the symmetric case with identical children. 

The severity of these two combinatorial problems depends on the event topology, as well 
as the type of signature objects. For example, there are cases where the first combinatorial 
problem is easily resolved, or even absent altogether. Consider the event topology of Fig. 3(a) 
with a lepton as the SM particle on each side. In this case, the partition is unique, and the 
upstream objects are jets, which can be easily identified [97]. Now consider the event topolo- 
gies of Figs. 3(b) and 3(c), with two opposite sign, same flavor leptons on each side. Such 
events result from inclusive pair production of heavier neutralinos in supersymmetry. By 
selecting events with different lepton flavors: e + e _ / u + / u~, we can overcome the first combina- 
torial problem above and uniquely associate the e + e~ pair with one cascade and the ^ \i~ 
pair with the other. However, the second combinatorial problem remains, as we still have 
to decide which of the two lepton pairs to associate with A = a and which to associate with 
A = b. Recall that the labels A = a and A = b are already attached to the child particles, 
which are distinguishable in the asymmetric case. In this paper we use the convention that 
A = a is attached to the lighter child particle: 



which also ensures that the M__ parameter defined in (3.9) is real. 

We can put this discussion in more formal terms as follows. The correct association of 
the visible particles with the corresponding children will yield 



M c (a) < M c (b) 



(3.26) 



M T 2(p T (a) ,P T (b) ; 



vrn, 



ia) ,m (b yM a) M b) ), 



(3.27) 



while the other, wrong association will give simply 



M T 2(p T (a) ,P T (b) ; 



vrn, 



ia) ,m (b yM b) M a) ). 



(3.28) 
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Figure 4: Unit-normalized Mt2 distributions for the event topology of Fig. 3(b). The mass spectrum 
is chosen as M c (a) = 100 GeV, M c (b) = 200 GeV and M p = 600 GeV. The test children masses are 
taken to be the true masses: = and Ms — M^ 1 . The dotted black distribution is the true 

Mt2 distribution, ignoring the combinatorial problem. The red histogram shows the distribution of the 
Mj^ variable defined in (3.29) while the blue histogram shows the distribution of the variable 
defined in (3.31). 

Both of these two M?2 values can be computed from the data, but a priori we do not know 
which one corresponds to the correct association. The solution to this problem is however 
already known [64,85]: one can conservatively use the smaller of the two 

M$ = min { M T2 (p T (a) ,p T w ;m (a) ,m (b) -M a) M b) ), M T2 (p t w ;m (a) ,m (b) ;M C (6) M a) ) } (3.29) 

in order to preserve the location of the upper Mt2 endpoint. This is illustrated in Fig. 4, 
where we show results for the event topology of Fig. 3(b) with a mass spectrum as follows: 
M c {a) = 100 GeV, M c (fc) = 200 GeV and M p = 600 GeV. The test children masses are taken 
to be the true masses: = and = Mc b \ The dotted black distribution is 

the unit-normalized true Mt2 distribution, where one ignores the combinatorial problem and 
uses the Monte Carlo information to make the correct association. The red histogram shows 
the unit-normalized distribution of the variable defined in (3.29). We see that the 

definition (3.29) preserves the corresponding endpoint: 

M Tf(ma x) = M T2( maX ). (3.30) 

Of course, we can also consider the alternative combination 

M$ = max{M T2 (p T Wp T ^m (a) ,^ (3.31) 
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whose unit-normalized distribution is shown in Fig. 4 with the blue histogram. One can 
see that some of the wrong combination entries in the histogram violate the original 

endpoint MT2(max) > yet there is still a well defined Mj^ endpoint 

M T2 { max)> M T2 { max)- (3.32) 

Strictly speaking, in our analysis in the next sections, we only need to study the 
endpoint (3.30), which contains the relevant information about the physical Mt2 endpoint. At 
the same time, with our convention (3.26) for the children masses, we only need to concentrate 
on the upper half > of the (Mc a \ M^) plane. However, for completeness we shall 
also present results for the endpoint (3.32), and we shall use the lower (M c (6) < M c (a) ) 

half of the (mJ°\ mJ 6 ^) plane to show those. Thus the Mt2 endpoint shown in our plots 
below should be interpreted as follows 

f M$, „ if a) <0 b) 



4. The simplest event topology: one SM particle on each side 

In this section, we consider the simplest topology with a single visible particle on each side of 
the event. We already introduced this example in Section 3.4, along with its event topology 
in Fig. 3(a). In Section 4.1 below we first discuss an asymmetric case with different children. 
Later in Section 4.2 we consider a symmetric situation with identical children masses. The 
mass spectra for these two study points are listed in Table 1. 



Spectrum 


Case 


M c (a) 




M p 


I 


Different children 


250 


500 


600 


II 


Identical children 


100 


100 


300 



Table 1: Mass spectra for the two examples studied in Sections 4.1 and 4.2. All masses are given in 
GcV. 



4.1 Asymmetric case 

Before we present our numerical results, it will be useful to derive an analytical expression 
for the asymmetric Mt2 endpoint (3.15) in terms of the corresponding physical spectrum 
of Table 1 and the two test children masses and Our result will generalize the 

corresponding formula derived in [68] for the symmetric case of M c (a) = M c (6) = M c and no 
upstream momentum {Pjjtm = 0). For the event topology of Fig. 3(a) the M.T2 endpoint is 
always obtained from the balanced solution and is given by [68] 

M T 2(max) (Mc, PuTM = 0) = /i ppc + ^ 2 ppc + . (4.1) 
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Here we made use of the convenient shorthand notation introduced in [85] for the relevant 
combination of physical masses 




AVc=^ 1- ^ • (4-2) 



The fi parameter defined in (4.2) is simply the transverse momentum of the (massless) visible 
particle in those events which give the maximum value of Mt2 [97]. Squaring (4.1), we can 
equivalently rewrite it as 



(M c , Putm = 0) = 2 4 pc + Ml + ^4 m| pc (m| P c + Ml) . (4.3) 

Now let us derive the analogous expressions for the asymmetric case M c (a) / M c (b) . Just 
like the symmetric case, the asymmetric endpoint M T2 ( ma x) also comes from a balanced 
solution and is given by 



M T2{max) (M c (a) , M c (6) , Putm = 0) = 2^ pc + M\ + y 4 p? ppc {j% pc + M\ ) + Mf , (4.4) 

where the parameters M 2 . and M 2 were already defined in (3.8) and (3.9), while p, ppc is now 
the geometric average of the corresponding individual fi ppc parameters 

_ (M p 2 -(M c (a) ) 2 )(M p 2 -(M c (fe) ) 2 ) 



-2 _ 

It is easy to check that in the symmetric limit 



MW M(a) ^ Uppc ^ m + Mc ^ M _ ^ 0> (4. 6 ) 

eq. (4.4) reduces to its symmetric counterpart (4.3), as it should. 

We are now ready to present our numerical results for the event topology of Fig. 3(a). 
We first take the asymmetric mass spectrum I from Table 1 and consider the case with no up- 
stream momentum, when formula (4.4) applies. Fig. 5 shows the corresponding Mt2 endpoint 
as a function of the two test children masses and . In panel (a) we present a three 
dimensional view, while in panel (b) we show a contour plot projection on the {M^\m^) 
plane (red contour lines) . On either panel, the green dot marks the true values of the children 
masses, Mc^ and Panel (b) also shows a gradient plot, where longer (shorter) arrows 

imply steeper (gentler) slope. The symmetric endpoint M T2 ( max )(M c , Putm = 0) of eq. (4.1) 
can be obtained by going along the diagonal orange line = in Fig. 5(b). We remind 
the reader that the endpoint Mj>2(max) plotted in Fig. 5 should be interpreted as in eq. (3.33). 

Fig. 5 illustrates the first basic property of the asymmetric Mt2 variable, which was 
discussed in Sec. 3.3.1. The M?2 endpoint allows us to find one relation between the two 
children masses and and the parent mass M p = M T2 (max) > an d in order to do so, we 
do not have to assume equality of the children masses, as is always done in the literature. The 
crucial advantage of our approach, in which we allow the two children masses to be arbitrary, is 
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Figure 5: Mx2(max) as a function of the two test children masses, Mc and M c , for the event 
topology of Fig. 3(a) with no upstream momentum (Pjjtm — 0), and the asymmetric mass spectrum 
I from Table 1: (M c (q) , M c (b) , M p ) = (250,500,600) GeV. We show (a) a three dimensional view and 
(b) contour plot projection on the (Mc a \ Mc) plane (red contour lines). The green dot marks the 
true values of the children masses. Panel (b) also shows a gradient plot, where longer (shorter) 
arrows imply steeper (gentler) slope. A kink structure is absent in this case. The symmetric endpoint 
MT2(max)(M c ) of eq. (4-1) can be obtained by going along the diagonal orange line Mc = M^ . 

its generality and model-independence. It allows us to extract the basic information contained 
in the Mt2 endpoint, without muddling it up with additional theoretical (and unproven) 
assumptions. 

Unfortunately, to go any further and determine each individual mass, we must make use 
of the additional properties discussed in Sees. 3.3.2 and 3.3.3. In the case of the simplest 
event topology of Fig. 3(a) considered here, they both require the presence of some upstream 
momentum [67,85]. As a proof of concept, we now reconsider the same type of events, but 
with a fixed upstream momentum of Pjjtm = 1 TeV. (The upstream momentum may be due 
to initial state radiation, or decays of heavier particles upstream.) The corresponding results 
are shown in Fig. 6. 

Fig. 6 demonstrates the second basic property of the asymmetric Mt2 variable discussed 
in Sec. 3.3.2. Unlike the result shown in Fig. 5(a), which was perfectly smooth, this time the 
M T2 (max) function in Fig. 6(a) shows a ridge, corresponding to the slope discontinuity marked 
with the black solid line in Fig. 6(b). The most important feature of the ridge is the fact 
that it passes through the green dot marking the true values of the children masses. Notice 
that applying the traditional symmetric Mt2 approach in this case will give a completely 
wrong result. If we were to assume equal children masses from the very beginning, we will be 
constrained to the diagonal orange line in Fig. 6(b). The Mt2 endpoint will then still exhibit 
a kink, but the kink will be in the wrong location. In the example shown in Fig. 6(b), we 
will underestimate the parent mass, while for the child mass we will find a value which is 
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Figure 6: The same as in Fig. 5 but with fixed upstream momentum of Put M = 1 TeV. The ridge 
structure (shown as the black solid line) is revealed by the sudden increase in the slope (gradient) in 
panel (b). Notice that the ridge goes through the true values of the children masses marked by the 
green dot. 

somewhere in between the two true masses M c {a) and M c (b) . 

Using the ridge information, we now know an additional relation among the children 
masses, which allows us to express all three masses in terms of a single unknown parameter 
9, as illustrated in Fig. 7(a). Let us choose to parametrize the ridge by the polar angle in the 



Using the ridge information from Fig. 6, we can then find all three masses as a function of 
9. The result is shown in Fig. 7(a). The mass of the lighter child is plotted in red, the 
mass Mc of the heavier child is plotted in blue, while the parent mass M p is plotted in black. 
With our convention (3.26) for the children masses, only values of 9 > 45° are physical, and 
the corresponding masses are shown with solid lines. The dotted lines in Fig. 7(a) show the 
extrapolation into the unphysical region 9 < 45°. 

Fig. 7(a) has some important and far reaching implications. For example, one may now 
start asking the question: Are there really any massive invisible particles in those events, 
or is the missing energy simply due to neutrino production [21]? The ridge results shown 
in Fig. 7(a) begin to provide the answer to that quite fundamental question. According to 
Fig. 7(a), for any value of the (still unknown) parameter 9, the two children particles cannot 
be simultaneously massless. This means that the missing energy cannot be simply due to 
neutrinos, i.e. there is at least one new, massive invisible particle produced in the missing 
energy events. At this point, we cannot be certain that this is a dark matter particle, but 
establishing the production of a WIMP candidate at a collider is by itself a tremendously 
important result. Notice that while we cannot be sure about the masses of the children, the 



(M c (a) ,M c (b) ) pi; 



ane: 




(4.7) 




Figure 7: (a) Particle masses obtained along the M T 2(max) ridge seen in Fig. 6. The ridge is 
parametrized by the angle 9 defined in (4--^)- The two children masses M { c a) {9) (in red) and M c (b) (6>) 
(in blue) as well as the parent mass M p (in black) are then plotted as a function of 9. In our con- 
vention (3.26) only values of 9 > 45° are physical, and the corresponding masses are shown with 
solid lines. Dotted lines show the extrapolation for 9 < 45°. (b) Contour plot of the quantity 
AM T2 ( max) {M^ a) , , P UTM = 1 TeV) defined in eq. (3.21), in the (Afi a) , M c (b) ) plane. This plot 
is obtained simply by taking the difference between Fig. 6(a) and Fig. 5(a). The solid black curve 
indicates the location of the Mx2(max) ridge. Only the point corresponding to the true children masses 
(the green dot) satisfies the Pjjtm invariance condition AM T 2( max ) = from eq. (3.23). 

parent mass M p is determined with a very good precision from Fig. 7(a): the function M p (9) 
is almost flat and rather insensitive to the particular value of 9 14 . 

Once we have proved that some kind of WIMP production is going on, the next imme- 
diate question is: how many such WIMP particles are present in the data - one or two? 
Unfortunately, the ridge analysis of Fig. 7(a) alone cannot provide the answer to this ques- 
tion, since the value of 9 is still undetermined. If 9 = 90°, one of the missing particles is 
massless, which is consistent with a SM neutrino. Therefore, if 6 were indeed 90°, the most 
plausible explanation of this scenario would be that only one of the missing particles is a 
genuine WIMP, while the other is a SM neutrino. On the other hand, almost any other value 
of 9 < 90° would guarantee that there are two WIMP candidates in each event. In that case, 
the next immediate question is: are they the same or are they different? Fortunately, our 
asymmetric approach will allow answering this question in a model-independent way. If is 
determined to be 45°, the two WIMP particles are the same, i.e. we are producing a single 
species of dark matter. On the other hand, if 45° < 9 < 90°, then we can be certain that 
there are not one, but two different WIMP particles being produced. 

14 Interestingly, for the example in Fig. 7(a), the maximum value of M p (6) happens to give the true parent 
mass M p , but we have checked that this is a coincidence and does not hold in general for other examples which 
we have studied. 
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We see that in order to completely understand the physics behind the missing energy sig- 
nal, we must determine the value of 6, i.e. we must find the exact location of the true children 
masses along the ridge. One of our main results in this paper is that this can be done by using 
the third Mt2 property discussed in Sec. 3.3.3. The idea is illustrated in Fig. 7(b), where we 
show a contour plot in the (mJ°\ M<^) plane of the quantity AM T2 ( max )(Mc a \ Putm) 
defined in eq. (3.21), for a fixed Putm = 1 TeV. This plot is obtained simply by taking the 
difference between Fig. 6(a) and Fig. 5(a). (A more practical method for obtaining this infor- 
mation was proposed in [100].) Recall that the function AMr 2 ( mm ) was introduced in order 
to quantify the Pjjtm invariance of the Mt2 endpoint, and it is expected that AM T2 ( m ax) 
vanishes at the correct values of the children masses (see eq. (3.23)). This expectation is 
confirmed in Fig. 7(b), where we find the minimum (zero) of the AM T2 ( m ax) function exactly 
at the right spot (marked with the green dot) along the M T2 ^ max ^ ridge. Thus the AM T2 ( mai ) 
function in Fig. 7(b) completely pins down the spectrum, and in this case would reveal the 
presence of two different WIMP particles, with unequal masses M c (a) / M c (fe) . Our analysis 
thus shows that colliders can not only produce a WIMP dark matter candidate and measure 
its mass, as discussed in the existing literature, but they can do a much more elaborate 
dark matter particle spectroscopy, as advertized in the title. In particular, they can probe 
the number and type of missing particles, including particles from subdominant dark matter 
species, which are otherwise unlikely to be discovered experimentally in the usual dark matter 
searches. 

4.2 Symmetric case 

While in our approach the two children masses M c (a) and M c (6) are treated as independent 
inputs, this, of course, does not mean that the approach is only valid in cases when the children 
masses are different to begin with. The techniques discussed in the previous subsection 
remain applicable also in the more conventional case when the children are identical, i.e. when 
colliders produce a single dark matter component. In order to illustrate how our method works 
in that case, we shall now work out an example with equal children masses. We still consider 
the simplest event topology of Fig. 3(a), but with the symmetric mass spectrum II from 
Table 1. We then repeat the analysis done in Figs. 5, 6, and 7 and show the corresponding 
results in Figs. 8, 9 and 10. 

The conclusions from this exercise are very similar to what we found earlier in Sec. 4.1 
for the asymmetric case. The Mt2 endpoint still provides one relation among the two chil- 
dren masses Mc^ and and the parent mass M p = M T2 r max \. This relation is shown in 
Fig. 8 (Fig. 9) for the case without (with) upstream momentum Putm- As seen in Fig. 8, 
in the absence of any upstream Putm, the function M p (M^ a \ M^) is smooth and reveals 
nothing about the children masses. However, the presence of upstream momentum signifi- 
cantly changes the picture and the function M p {M^ a \ M^) again develops a ridge, which is 
clearly visible 15 in both the three-dimensional view of Fig. 9(a), as well as the gradient plot 
in Fig. 9(b). The ridge information now further constrains the children masses to the black 

15 We caution the reader that here we are presenting only a proof of concept. In the actual analysis the ridge 
may be rather difficult to see, for a variety of reasons - detector resolution, finite statistics, combinatorial and 
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Figure 8: The same as Fig. 5, but for the symmetric mass spectrum II from Table 1, i.e. 
(M c (a) ,M c (b) ,M p ) = (100,100,300) GeV. 




Figure 9: The same as Fig. 6 but for the symmetric mass spectrum II from Table 1, i.e. 
(Afi a) , M c (b) , M p ) = (100,100,300) GeV. 

solid line in Fig. 9(b), leaving only one unknown degree of freedom. Parametrizing it with the 
polar angle 9 as in (4.7), we obtain the spectrum as a function of 9, as shown in Fig. 10(a). 
Once again we find the fortuitous result that in spite of the remaining arbitrariness in the 
value of 9, the parent mass M p is very well determined, since M p (9) is a very weakly varying 
function of 9. Furthermore, both Fig. 9(a) and Fig. 9(b) exhibit a high degree of symmetry 

SM backgrounds, etc. Nevertheless, we expect that the ridge will be just as easily observable as the traditional 
kink in the symmetric Mti endpoint. If the kink can be seen in the data, the ridge can be seen too, and there 
is no reason to make the assumption of equal children masses. Conversely, if the kink is too difficult to see, 
the ridge will remain hidden as well. 
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Figure 10: The same as in Fig. 7 but for the symmetric mass spectrum II from Table 1, i.e. 
(M c (a) ,M c (h) ,M p ) = (100,100,300) GeV. Notice that, in contrast to Fig. 7, the minimum of the 
AMT2(max) function is now obtained at Mc = Mc b \ implying that the two missing particles are 
the same. 

under M c (a) <-> M^ h \ which is a good hint that the children are in fact identical. This suspi- 
cion is confirmed in Fig. 10(b), where we find that the Putm dependence disappears at the 
symmetric point = Mc = 100 GeV, revealing the true masses of the two children. 

In the two examples considered so far in Sections 4.1 and 4.2, we used a fixed finite value 
of the upstream transverse momentum Putm = 1 TeV, which is probably rather extreme - 
in realistic models, one might expect typical values of Putm on the order of several hundred 
GeV. However, things begin to get much more interesting if one were to consider even larger 
values of Putm- On the one hand, the ridge feature becomes sharper and easier to observe 
[85]. More importantly, the ridge structure itself is modified, and a second set of ridgelines 
appears 16 at sufficiently large Putm- All ridgelines intersect precisely at the point marking 
the true values of the children masses, thus allowing the complete determination of the mass 
spectrum by the ridge method alone. This procedure was demonstrated explicitly in Ref. [96], 
which investigated the extreme case of Putm = oo for a study point with different parents 
and identical children. The assumption of Putm = oo justified the use of a "decoupling 
argument", in which the two branches A = a and A = b are treated independently, allowing 
the derivation of simple analytical expressions for the Mt2 endpoint [96]. In Appendix A 
we reproduce the analogous analytical results at Putm — > oo for the case of interest here 
(identical parents and different children) and study in detail the Putm dependence of the 
ridgelines. Unfortunately, we find that the values of Putm necessary to reveal the additional 
ridge structure, are too large to be of any interest experimentally. On the positive side, the 
Putm invariance method discussed in Sec. 2.3.3 does not require such extremely large values 
of Putm and can in principle be tested in more realistic experimental conditions. 

16 A keen observer may have already noticed a hint of those in Figs. 7(b) and 10(b). 
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4.3 Mixed case 



For simplicity, so far in our discussion we have been studying only one type of missing energy 
events at a time. In reality, the missing energy sample may contain several different types 
of events, and the corresponding Mt2 measurements will first need to be disentangled from 
each other. 

For concreteness, consider the inclusive pair production of some parent particle Xpi which 
can decay either to a child particle Xa of mass M c (a) , or a different child particle Xb of mass 
Mc b \ Let the corresponding branching fractions be B a and B^, i.e. B a = B{x p — > Xa) and 
Bb = B(x p — > Xb)- Furthermore, let Xb decay invisibly 17 to Xa- Such a situation can be 
easily realized in supersymmetry, for example, with the parent being a squark, a slepton, or 
a gluino, the heavier child Xb being a Wino-like neutralino x% and the lighter child Xa being a 
Bino-like neutralino x.i- The heavier neutralino has a large invisible decay mode x\ ~^ x\ v ^i 
if its mass happens to fall between the sneutrino mass and the left-handed slepton mass: 

m,<m^<m 1l . 

Let us start with a certain total number of events N pp in which two parent particles Xp 
have been produced. Then the missing energy sample will contain iV&& = N pp B^ symmetric 
events where the two children are Xb and Xb, N aa = N pp B^ symmetric events where the two 
children are Xa and Xa, and N a b = 2N pp B a Bi J asymmetric events where the two children are 
Xa and Xb- How can one analyze such a mixed event sample with a single Mt2 variable? 

The black histogram in Fig. 11 shows the unit-normalized Mt2 distribution for the whole 
(mixed) event sample (for convenience, we do not show the zero bin [100]). For this plot, we 
used the asymmetric mass spectrum I from table 1: = 250 GeV, = 500 GeV and 



M p = 600 GeV, and chose zero test masses for the children = M c u> = 0. For definiteness, 
we fixed equal branching fractions B a = B = 50%, so that the relative normalization of the 
three individual samples is N aa : : N a jj = 1:1:2. Fig. 11 shows that the observable 
Mt2 distribution is simply a superposition of the Mt2 distributions of the three individual 
samples XaXa, XaXb and XbXb, which are shown with the red, blue and green histograms, 
correspondingly. Each individual sample exhibits its own Mt2 endpoint, marked with a 
vertical arrow, which can also be seen in the combined Mt2 distribution. Using eq. (4.4), the 
three endpoints are found to be 



XaXa 



M, 



(aa) 



T2(max) 



(0,0,0) = M p 



M, 



(aY 



496 GeV, 



(4.8) 



XaXb ^ M^l^OAO) = M p 



\ 



M, 



(a)' 



M. 



0>y 



M n 



301 GeV, (4.9) 



XbXb ^ M^} max) (0,0,0) = M p 



Mc 



or 



M n 



183 GeV. 



(4.10) 



7 If Xb decays visibly, then the respective types of events can in principle be sorted by their signature. 
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M T2 (GeV) 

Figure 11: Unit-normalized, zero-bin subtracted Mt2 distribution (black histogram) for the full mixed 
event sample, as well as the individual components XaXa (red), XaXb (blue) and XbXb (green). We 
took zero test masses for the children Mc°^ = M^ 1 = and equal branching fraction for the parents 
B a — Bb — 50%. The mass spectrum is taken from the asymmetric study point I in Table 1 with 
M c (q) = 250 GeV, M c (6) = 500 GeV and M p = 600 GeV. The three arrows indicate the expected 
endpoints for each individual component in the sample. 

Now suppose that all three endpoints (4.8-4.10) are seen in the data. Their interpretation 
is far from obvious, and in fact, there will be different competing explanations. If one insists 
on the single missing particle hypothesis, there can be only one type of child particle, and the 
only way to get three different endpoints in Fig. 11 is to have production of three different 
pairs of parent particles, each of which decays in exactly the same way. Since the three parent 
masses are a priori unrelated, one does not expect any particular correlation among the three 
observed endpoints (4.8-4.10). Now consider an alternative explanation where we produce a 
single type of parents, but have two different children types. This situation also gives rise 
to three different event topologies, with three different Mt2 endpoints, as we just discussed. 
However, now there is a predicted relation among the three Mt2 endpoints, which follows 
simply from eqs. (4.8-4.10): 

If the parents are the same and the children are different, this relation must be satisfied. 
If the parents are different and the children are the same, a priori there is no reason why 
eq. (4.11) should hold, and if it does, it must be by pure coincidence. The prediction (4.11) 
therefore is a direct test of the number of children particles. Another test can be performed if 
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we could estimate the individual event counts N aa , N ab and N^, although this appears rather 
difficult, due to the unknown shape of the Mt2 distributions in Fig. 11. In the asymmetric 
example discussed here, we have another prediction, namely 

N ab = 2^/N aa N bb , (4.12) 

which is another test of the different children hypothesis. Notice that eq. (4.12) holds re- 
gardless of the branching fractions B a and B b , although if one of them dominates, the two 
endpoints which require the other (rare) decay may be too difficult to observe. 

Of course, the ultimate test of the single missing particle hypothesis is the behavior of the 
intermediate M?2 endpoint in Fig. 11 corresponding to the asymmetric events of type XaXb- 
Applying either one of the two mass determination methods discussed earlier in Figs. 7 and 
10, we should find that M^) > is a result of asymmetric events, indicating the simultaneous 
presence of two different invisible particles in the data. 

5. A more complex event topology: two SM particles on each side 

In this section, we consider two more examples: the off-shell event topology of Fig. 3(b) is 
discussed in Sec. 5.1, while the on-shell event topology of Fig. 3(c) is discussed in Sec. 5.2. (For 
simplicity, we do not consider any Putm in this section.) Now there are two visible particles 
in each leg, which form a composite visible particle of varying mass m^x). In general, by 
studying the invariant mass distribution of m^, one should be able to observe two different 
invariant mass endpoints, suggesting some type of an asymmetric scenario. 

5.1 Off-shell intermediate particle 

Here we concentrate on the example of Fig. 3(b). Since the intermediate particle is offshell, 
the maximum kinematically allowed value for mpg is given by eq. (3.24). 

Recall that for the simple topology of Fig. 3(a) discussed in the previous section, the Mt2 
endpoint (4.4) always corresponded to a balanced solution. More precisely, the Mt2 variable 
was maximized for a momentum configuration in which Mt2 was given by the balanced 
solution (3.7). However, in this section we shall find that for the more complex topologies of 
Figs. 3(b) and 3(c), the Mt2 endpoint may result from one of four different cases altogether: 
two different balanced solutions, which we shall label as B and B', or the unbalanced solutions 
Ua and Ub discussed in Sec. 3.2. Depending on the type of solution giving the endpoint 
M T2 ( max ), the (mJ°^ , ) parameter plane divides into the three regions 18 shown in Fig. 12. 
The green dot in Fig. 12 denotes the true children masses in this parameter space. Within each 
region, we show the relevant momentum configuration for the visible particles (red arrows) 
and the children particles (blue arrows) in each leg (a or b). The momenta are quoted in the 
"back-to-back boosted" (BB) frame [68], in which the two parents are at rest. The length 
of an arrow is indicative of the magnitude of the momentum. A blue dot implies that the 

18 The fourth case of the B' balanced solution happens to coincide with the two unbalanced solutions along 
the boundary between Ua and Ub. 
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Figure 12: The four regions in the (_Mi° , Mc ) parameter plane leading to the four different types 
of solutions for the Mt2 endpoint, for the off-shell event topology of Fig. 3(b). The green dot marks 
the true location of the two children masses. Within each region, we indicate the relevant momentum 
configuration for the visible particles (red arrows) and the children particles (blue arrows) in each leg 
(a or b). The momenta are quoted in the "back-to-back boosted" (BB) frame [68], in which the two 
parents are at rest. A blue dot implies that the corresponding daughter is at rest and therefore the two 
visible particles are emitted back-to-back. The two balanced solutions are denoted as B and B' , while 
the two unbalanced solutions are Ua and Ub. The black solid lines represent phase changes between 
different solution types and delineate the expected locations of the ridges in the M T 2(max) function 
shown in Fig. 13. 



corresponding daughter is at rest and therefore the two visible particles are emitted back- 
to-back. The two balanced solutions are denoted as B and B', while the two unbalanced 
solutions are Ua and Ub. The black solid lines represent phase changes between different 
solution types and delineate the expected locations of the ridges in the M^imax) function 
shown in Fig. 13 below. Perhaps the most striking feature of Fig. 12 is that the three (in fact, 
all four) regions come together precisely at the green dot marking the true values of the two 
children masses. The boundaries of the regions shown in Fig. 12 will manifest themselves as 
the locations of the ridges (i.e. gradient discontinuities) in the M T 2(max) function. Therefore, 
we expect that by studying the ridge structure and finding its "triple" point, one will be able 
to completely determine the mass spectrum. 

We shall now give analytical formulas for the Mt2 endpoint in each of the four regions 
of Fig. 12. We begin with the two balanced solutions B and B' , for which the event-by-event 
balanced solution for Mt2 is given by eq. (3.7). In the parameter space region of Fig. 12 
which is adjacent to the origin, we find the balanced configuration B, in which all visible 
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particles have the same direction in the BB frame. As a result, we have 

m (a) = m (b) = (5.1) 

and 

= 2M| • (5 ' 2) 

Substituting eqs. (5.1) and (5.2) in the balanced Mti solution (3.7), where we should take 
the plus sign, we obtain 



M? 2(max) (M^\M^)] = 2/4 c + Ml + ^4 4 pc (flj pc + Ml) + Mi, (5.3) 

which we recognize as the balanced solution (4.4) found for the decay topology of Fig. 3(a). 

Moving away from the origin in Fig. 12, we find a second balanced solution B' along 
the boundary of the unbalanced regions Ua and Ub. In this case the visible particles are 
back-to-back, and their invariant mass is maximized: 

m {x) =M P - M C W , (5.4) 

and correspondingly 



A, 



M p - M c (°>) (M p - AfW) . (5.5) 

Substituting eqs. (5.4) and (5.5) in the balanced Mt2 solution (3.7), we obtain the .B'-type 
Mt2 endpoint as 

M* 2{max) (M^\m} 2 = ( Mp -MP) (m p -M<«) +Ml+ 2Mp M ^ M tf ] Ml- 

(5.6) 

The corresponding formulas for the unbalanced cases Ua and Ub are obtained by taking 
the maximum value for the invariant mass of the visible particles in the corresponding decay 
chain: 

m (a) = mV^ x = M p - M c (a ) for region (Ua) , (5.7) 
m{b) = m ™ x = M p - M^ for region (Ub) . (5.8) 

The corresponding formula for M T2 ( m ax) 1S then given by 

< 2 W) )=M P - M(») + MW , (5.9) 
M^ {max) (M^) =M P - M( fc ) + M^ . (5.10) 

One can now use the analytical results (5.3), (5.6), (5.9) and (5.10) to understand the 
ridge structure shown in Fig. 12. For example, the boundary between the B and Ua regions 
is parametrically given by the condition 

M* 2{max) (M^,M^) = M^ (max) (M^) , (5.11) 
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Figure 13: The same as Fig. 5, but for the off-shell event topology of Fig. 3(b). We use the mass 
spectrum from the example in Fig. 4-' Mc — 100 GeV, Mc = 200 GeV and Al p = 600 GeV and for 
simplicity consider only events with Putm = 0. 

while the boundary between the B and Ub regions is parametrically given by 

M^ {max) (Mt\W) = M™ {max) (MW) . (5.12) 
On the other hand, the boundary 

<2W*i 0) ) = M™ {max) (MW) (5.13) 

between the two unbalanced regions Ua and Ub is quite interesting. The parametric equation 
(5.13) is nothing but a straight line in the {M^\Mc) plane: 

M c ( fe ) = M C W - + AfW , (5.14) 

as seen in Fig. 12. 

It is now easy to understand the triple point structure in Fig. 12. The triple point is 
obtained by the merging of all three boundaries (5.11), (5.12) and (5.13), i.e. when 

M^ 2{max) (Mt\M^) = M T B ; M (MW,Mf)) = M^ mffl) (MW) = M^ M #) . 

(5-15) 

It is easy to check that = M^ 1 and Mj = identically satisfy these equations, 

thereby proving that the triple intersection of the boundaries seen in Fig. 12 indeed takes 
place at the true values of the children masses. 

These results are confirmed in our numerical simulations. In Fig. 13 we present (a) a 
three dimensional view and (b) a gradient plot of the ridge structure found in events with the 
off-shell topology of Fig. 3(b). The mass spectrum for this study point was fixed as in Fig. 4, 
namely M c (a) = 100 GeV, M c (fe) = 200 GeV and M p = 600 GeV. Since the ridge structure for 
this topology does not require the presence of upstream momentum, for simplicity we consider 
only events with Pjjtm = 0. The ridge pattern is clearly evident in Fig. 13(a), which shows 
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a three-dimensional view of the Mt2 endpoint function M T2 ( m ax)(Mc , Mc ). It is even 
more apparent in Fig. 13(b), where one can see a sharp gradient change along the ridge lines: 
in regions Ua and Ub, the corresponding gradient vectors point in trivial directions (either 
horizontally or vertically), in accord with eqs. (5.9)-(5.10). On the other hand, the gradient 
in region B is very small, and the Mt2 endpoint function is rather flat. The green dot marks 
the location of the true children masses (M c (a) = 100 GeV, M c (6) = 200 GeV) and is indeed 
the intersection point of the three ridgelines. As expected, the corresponding M T 2( m ax) at 
that point is the true parent particle mass M p = 600 GeV. 

At this point, it is interesting to ask the question, what would be the outcome of this 
exercise if one were to make the usual assumption of identical children, and apply the tradi- 
tional symmetric Mt2 to this situation. The answer can be deduced from Fig. 13(b), where 
the diagonal orange dotdashed line corresponds to the usual assumption of = In 
that case, one still finds a kink, but at the wrong location: in Fig. 13(b) the intersection of the 
diagonal orange line and the solid black ridgeline occurs at = Mc^ = 65.3 GeV and the 
corresponding parent mass is M p = 565.3 GeV. Therefore, the traditional kink method can 
easily lead to a wrong mass measurement. Then the only way to know that there was some- 
thing wrong with the measurement would be to study the effect of the upstream momentum 
and see that the observed kink is not invariant under Putm- 

We should note that, depending on the actual mass spectrum, the two-dimensional ridge 
pattern seen in Figs. 12 and 13(b) may look very differently. For example, the balanced 
region B may or may not include the origin. One can show that if 

Mp < ^ ( Mc(6) + M M ° a) ) 2 + ( M ° b) ) 2 ) » ( 5 - 16 ) 

the boundary between B and Ua does not cross the ax is- In this case the diagonal line 
in Fig. 13(b) does not cross any ridgelines and the traditional Mt2 approach will not produce 
any kink structure, in contradiction with one's expectations. This exercise teaches us that 
the failsafe approach to measuring the masses in missing energy events is to apply from the 
very beginning the asymmetric Mt2 concept advertized in this paper. 

5.2 On-shell intermediate particle 

Our final example is the on-shell event topology illustrated in Fig. 3(c). Now there is an 
additional parameter which enters the game — the mass of the intermediate particle 

in the A-th decay chain. As a result, the allowed range of invariant masses for the visible 
particle pair on each side is limited from above by eq. (3.25). 

In this case we find that the Mt2 endpoint exhibits a similar phase structure as the one 
shown in Fig. 12. One particular pattern is illustrated in Fig. 14, which exhibits the same four 
regions B, B', Ua and Ub seen in Fig. 12. The difference now is that region B' is considerably 
expanded, and as a result, region B does not have a common border with regions Ua and 
Ub any more. The triple point of Fig. 12 has now disappeared and the correct values of the 
children masses now lie somewhere on the border between regions B and B', but their exact 
location along this ridgeline is at this point unknown. 
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(b): - 




M 



(a) 



Figure 14: The same as Fig. 12 but for the onshell scenario illustrated in Fig. 3(c). 

Just like we did for the off-shell case in Sec. 5.1, we shall now present analytical formulas 
for the Mt2 endpoint in each region of Fig. 14. In the balanced region B, we find the same 
results (5.1-5.3) as in the off-shell case considered in the previous Section 5.1. The other 
balanced region B' is characterized by 



m 



(A) 



m 



(A) 



(5.17) 



where rri^ x is given by eq. (3.25), and 



An 



M 2 



+ 




(5.18) 



The formula for the endpoint ^^(max) * n re gi° n B' is then simply obtained by substituting 
(5.17) and (5.18) into the balanced solution (3.7). 

Finally, the Mt2 endpoint in the unbalanced regions U a and Ub is given by 



M T2(max) 
lvl T2(max) 



(Mi 6 )) 



m 



in 



(a) 

ma,', 
(b) 



+ M c (a) , 



(5.19) 
(5.20) 



where m™ x and m^ x are given by eq. (3.25). 

In Fig. 15 we present our numerical results in this on-shell scenario. The mass spectrum 



is fixed as: M, 



(a) 



100 GeV, M, 



CO 



200 GeV, M 



(«) 



M 



(&) 



550 GeV and M p = 1 TeV, 
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Figure 15: The same as in Fig. 13 but for the onshell scenario of Fig. 3(c), with a mass spectrum 
M c (q) = 100 GeV, M c (b) = 200 GeV, M^ a) = = 550 GeV and M p = 1 TeV. 

and we still do not include the effects of any upstream momentum. Fig. 15(a) shows the 
three-dimensional view of the Mra endpoint function M^2(max){^c a \ Mi ), which exhibits 
three different sets of ridges, which are more easily seen in the gradient plot of Fig. 15(b). 
As usual, the green dot marks the true children masses. Fig. 15(b) shows that the ridgeline 
separating the two balanced regions B and B' does go through the green dot and thus reveals 
a relationship between the two children masses, leaving the ridgeline parameter 9 as the only 
remaining unknown degree of freedom. However, unlike the off-shell case of Sec. 5.1, now 
there is no special point on this ridgeline, and we cannot completely pin down the masses by 
the ridge method. Thus, in order to determine all masses in the problem, one must use an 
additional piece of information, for example the visible invariant mass endpoint (3.25) or the 
Putm invariance method suggested in Sec. 3.3.3. 

6. Summary and conclusions 

Cosmological observations hint towards the existence of one or more hypothetical dark matter 
particles. The start of the Large Hadron Collider may offer an unique opportunity to produce 
and study dark matter in a high-energy experimental laboratory. Unfortunately, the dark 
matter signatures at colliders always involve missing transverse energy. Such events will 
be quite challenging to fully reconstruct and/or interpret. All previous studies have made 
(either explicitly or implicitly) the assumption that each event has two identical missing 
particles. Our main point in this paper is that this assumption is unnecessary, and by suitable 
modifications of the existing analysis techniques one can in principle test both the number 
and the type of missing particles in the data. Our proposal here was to modify the Cambridge 
Mt2 variable [50] by treating each children mass as an independent input parameter. In this 
approach, one obtains the M?2 endpoint M^lmax) as a function of the two children masses 
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Mc and Mc , and proceeds to study its properties. The two most important features of the 
thus obtained function M T2 (max)(Mc a \ M^), identified in this paper, were the following: 

• The function M T2 (max){Mc a \ M^) exhibits a ridge structure (i.e. a gradient discon- 
tinuity), as illustrated with specific examples in Figs. 6, 9, 13 and 15. The point 
corresponding to the correct children masses always lies on a ridgeline, thus the ridge- 
lines provide a model-independent constraint among the children masses, just like the 
Mt2 endpoint provides a model-independent constraint on the masses of the child(ren) 
and the parent. 

• In general, the Mt2 endpoint function also depends on the value of the upstream 
transverse momentum in the event: Mrp 2 ^ max ^(Mc a \ Mc b \ Ptjtm)- However, the Pjjtm 
dependence disappears completely for precisely the right values of the children masses, 
as seen in the examples of Figs. 7(b) and 10(b). This provides a second, quite general 
and model-independent, method for measuring the individual particle masses in such 
missing energy events. 

Before we conclude, we shall discuss a few other possible applications of the asymmetric 
Mt2 idea, besides the examples already considered in the paper. 

1. Invisible decays of the next-to-lightest particle. Most new physics models introduce some 
new massive and neutral particle which plays the role of a dark matter candidate. Often 
the very same models also contain other, heavier particles, which for collider purposes 
behave just like a dark matter candidate: they decay invisibly and result in missing 
energy in the detector. For example, in supersymmetry one may find an invisibly 
decaying sneutrino vi — > veXv m UED one finds an invisibly decaying KK neutrino 
u\ —7- 1/71, etc. These scenarios can easily generate an asymmetric event topology. For 
example, consider the strong production of a squark (q) pair, as illustrated in Fig. 16(a). 
One of the squarks subsequently decays to the second lightest neutralino x%-, which in 
turn decays to the lightest neutralino Xi by emitting two SM fermions X2 ~~ ^ ^ + ^~Xi 
(or X2 ~~ ^ JJXi)- The other squark decays to a chargino xf > which then decays to a 
sneutrino as xf ~~ ^ ^^i- Since vi can only decay invisibly, we obtain the asymmetric 
event topology outlined with the blue box in Fig. 16(a). The two squarks are the 
parents, the lightest neutralino x\ is the first child, and the sneutrino i>i is the second 
child. 

2. Applying Mt2 to an asymmetric subsystem. One can also apply the Mt2 idea even to 
events in which there is only one (or even no) missing particles to begin with. Such 
an example is shown in Fig. 16(b), where we consider ti production in the dilepton or 
semi-leptonic channel. In the first leg we can take bt as our visible system and the 
neutrino vg as the invisible particle, while in the other leg we can treat the 6-jet as the 
visible system and the VF-boson as the child particle. In this case, there still should be 
a ridge structure revealing the true t, W and v masses. 
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Figure 16: Event topology for the two examples discussed in Section 6. The black solid lines represent 
SM particles which are visible in the detector while red solid lines represent particles at intermediate 
sages. The missing particles are denoted by dotted lines, (a) Squark pair production with decay chains 
terminating in two different invisible particles (\i and vi, correspondingly). In this case decays 
invisibly, (b) The subsystem Mt2 variable applied to ti events. The W -boson in the lower leg is treated 
as a child particle and can decay either hadronically or leptonically. 

3. Multi- component dark matter. Of course, the model may contain two (or more) different 
genuine dark matter particles [104-112], whose production in various combinations will 
inevitably lead at times to asymmetric event topologies. 

In conclusion, our work shows that the Mt2 concept can be easily generalized to decay 
chains terminating in two different daughter particles. Nevertheless, the methods discussed in 
this paper allow to extract all masses involved in the decays, at least as a matter of principle. 
We believe that such methods will prove extremely useful, if a missing energy signal of new 
physics is seen at the Tevatron or the LHC. 
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A. Appendix: The asymmetric Mti in the limit of infinite Putm 

In this appendix we revisit our previous two examples from Sections 4.1 and 4.2, this time 
considering the infinitely large Putm limit [96]. While this situation is impossible to achieve 
in a real experiment, its advantage is that it can be treated by analytical means. In the 
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Figure 17: The parameter plane of test children masses squared, divided into the four different regions 
IZi used to define the Mt2 endpoint function (A.l). Their common boundaries Lij are parametrically 
defined in eqs. (A.6-A.9). The black dot corresponds to the true values of the children masses. 



Pjjtm — > oo limit, the "decoupling argument" of Ref. [96] holds, and one finds the following 
analytical expression for the Mt2 endpoint as a function of the two test children masses M c (a) 
and M r (6) : 



T2(max) 



oo 



Ml - (M c (a) )2 + (M c (a) )2, if (M c w, M c w ) G n lt 



(«) ii>( 6 )> 



Ml - (M c (b) )2 + (M c (b) )2, if (M c (a) , M c {6) ) e TZ 2 , 



M, 



(b) 



M 



JEjMp, 



M, 



(a) 



if (M c (a) ,M c (6) ) G^ 3 , 
if (M c (a) ,M c W ) G^ 4 , 



(A.l) 



where the four defining regions TZi, {i = 1, ... ,4) are shown in Fig. 17 and are defined as 
follows: 



Til : M c (6) < \J (M c (6) ) 2 - (M c (a) ) 2 + (M c (a) ) 2 A M^ < , 
K 2 : \l (Mi b) ) 2 - (M c (a) ) 2 + (M c (a) ) 2 < M c ( fe ) < , 

^ 4 : M„ (a) < £ (a) A M„ (b) 




(A.2) 
(A.3) 

(A.4) 
(A.5) 
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Since the functional expression for Mx 2 ( m ax) within each region TZi is different, there is in 
general a gradient discontinuity when crossing from one region into the next. Therefore, the 
ridges on the M T2 ( max ) hypersurface will appear along the common boundaries of the four 
regions TZi. Let us denote by Lij the boundary between regions TZi and TZj. As indicated in 
Fig. 17, each is a straight line in the parameter space of the children test masses squared 
and is given by 

L 12 : (M c (b) ) 2 = (M^) 2 - (M^) 2 + (M^) 2 , AfW < ; (A.6) 
L 23 : AfW = M C W , M c (a) < ; (A.7) 
M {h) 

L 34 : = -4- M« , Afto > ; (A.8) 

Mr 



L u : Afto = , < . (A.9) 

As seen in Fig. 17, all four lines meet at the true children mass point = Mc a \ 

= Mc b \ where in turn the Mt2 endpoint M T2 (max) gives the true parent mass M p , in 
accordance with eq. (3.16). 

With those preliminaries, we are now in a position to revisit our two examples from 
Sections 4.1 and 4.2. Figs. 18 and 19 are the corresponding analogues of Figs. 6 and 9 in the 
case of infinite Putm- Comparing with our earlier results, we notice both quantitative and 
qualitative changes in the ridge structure. First, the smooth ridge in Fig. 6(b) (Fig. 9(b)) has 
now been deformed into two straight line segments, one horizontal (L23) and the other vertical 
(L14), which meet at an angle of 90° precisely at the true values of the children masses. More 
importantly, Figs. 18 and 19 now exhibit another pair of ridges L\ 2 and L34 (plotted in red in 
Figs. 18(b) and 19(b)), which were absent from the earlier figures in Section 4. The system 
of four ridges seen in Figs. 18(a) and 19(a) is very similar to the crease structure observed in 
Ref. [96] . We thus confirm the result of Ref . [96] that in the infinite Putm limit there exist 
four different ridges, whose common intersection point reveals the true masses of the parent 
and children particles. 

At this point it is instructive to contrast the two sets of ridgelines: L23 and L14 (shown 
in Figs. 18(b) and 19(b) in black) versus L\ 2 and L34 (shown in Figs. 18(b) and 19(b) in 
red). The boundaries L23 and L14 separate the union of regions TZ\ and TZ 2 from the union 
of regions TZ3 and TZ&. Along those boundaries, we observe a transition in the configuration 
of visible momenta which yields the maximum possible value of M?2- More precisely, in 
regions TZ\ and TZ 2 we find that the visible momenta for M T2 ^ max ) are parallel to the 
direction of the upstream momentum Putm, while in regions TZ3 and TZ4 we find that 
are anti-parallel to Pjjtm- This fact remains true even at finite values of Putm, which is 
why the ridgelines L23 and L14 could also be seen in the earlier plots from Sec. 4 at finite 
Putm = 1 TeV. 

On the other hand, the ridgelines L\ 2 and L 3 4 shown in red in Figs. 18(b) and 19(b) 
are due to the "decoupling argument" [96], which is strictly valid only in the infinite Putm 
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Figure 18: The same as Fig. 6 but for Putm — > oo. 




Figure 19: The same as Fig. 9 but for Putm — > oo. 

limit. This is why these ridges become apparent only at very large values of Putm, and are 
gradually smeared out at smaller Putm- 

The evolution of the ridge structure as a function of Putm is shown in Figs. 20 and 21. 
In order to compare the sharpness of the four ridges, we choose to vary the test children 
masses and m| along a circle centered on their true values and with a fixed radius R. 
Such a circle is guaranteed to cross all four ridges, and can be parameterized in terms of an 
angular coordinate 4> as follows 

M c ( a ) (<f>) = + R cos <p, (A. 10) 

M c ( fe ) (0) = M c (b) + R sin <\>. (A. 11) 
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Figure 20: A study of the sharpness of the Mt2 ridge for the example considered in Sec. 4-1- The 
event topology is that of Fig. 3(a) and the mass spectrum is = 250 GeV, Mc = 500 GeV and 

M p = 600 GeV. We plot the asymmetric Mt2 endpoint M T2 ( ma x){Mc a \4>), Mc b \(j)), Putm), as a 
function of the angular variable <f> parameterizing the circle of radius R defined in eqs. (A. 10, A. 11). 
The radius R of the circle is taken to be R — 50 GeV in panel (a) and R — 5 GeV in panel (b). We 
present results for four different choices of the upstream momentum Putm o>s labelled in the plot. 




Figure 21: The same as Fig. 20, but for the example considered in Sec. ^.2, where the input mass 
spectrum is fixed as M c (a) = 100 GeV, M c (b) = 100 GeV and M p = 300 GeV. 



Then in Fig. 20 (Fig. 21) we plot the asymmetric M T2 endpoint M T2 ( ma x) (Mc (0), Mc (<f>), Putm), 
as a function of the angular variable <f>, for the case of mass spectrum I studied in Section 4.1 
(mass spectrum II studied in Section 4.2). The radius R is taken to be R = 50 GeV in panels 
(a) and R = 5 GeV in panels (b). We present results for four different choices of the up- 
stream momentum: Putm = 100 GeV (black lines), Putm = 1 TeV (blue lines), Putm = 4 
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TeV (magenta lines), and Pjjtm = °o (red lines). Notice that the red lines at Pjjtm = oo 
in Figs. 20 and 21 are directly correlated to the three-dimensional plots of Figs. 18 and 19, 
while the blue lines at Pjjtm = 1 TeV in Figs. 20 and 21 are directly correlated to the 
three-dimensional plots of Figs. 6 and 9. 

Each one of the previously discussed ridges manifests itself as a kink in Figs. 20 and 
21. Indeed, the red lines for Pjjtm = oo reveal four clear kinks, which (from left to right) 
correspond to the ridgelines L34, L23, and L14. Using eqs. (A.6-A.9), it is easy to find 
the expected location of each kink in the Pjjtm -> 00 limit: 4> = {63.4°, 180°, 204.9°, 270°} 
for Fig. 20(a), <p = {63.4°, 180°, 206.4°, 270°} for Fig. 20(b), and <j> = {45°, 180°, 225°, 270°} 
for Figs. 21(a) and 21(b). However, as the upstream momentum is lowered to more realistic 
values, the kinks gradually wash out, albeit to a different degree. As anticipated from our 
earlier results, the smearing effect is quite severe for L 34 and L 12 , and by the time we reach 
Pjjtm = 1 TeV, those two kinks have completely disappeared. On the other hand, L23 and 
L14 are affected to a lesser degree by the smearing effect and are still visible at Pjjtm = 1 
TeV, but by Pjjtm = 100 GeV they are essentially gone as well. Notice that the variation 
in Pjjtm affects not only the sharpness of the kinks, but also their location. This was to be 
expected, since we already saw that the shape of the ridge is different at Pjjtm = 1 TeV and 
Pjjtm = co: compare the black ridge lines in Figs. 6(b) and 9(b) to those in Figs. 18(b) and 
19(b). Finally, as a curious fact we notice that the results shown in panels (a) and panels (b) 
of Figs. 20 and 21 are approximately related by a simple scaling with a constant factor. 
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